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We introduce exactly solvable models of interacting (Majorana) fermions in d > 3 spatial dimen¬ 
sions that realize a new kind of topological quantum order, building on a model presented in ref. 

[1]. These models have extensive topological ground-state degeneracy and a hierarchy of point-like, 
topological excitations that are only free to move within sub-manifolds of the lattice. In particular, 
one of our models has fundamental excitations that are completely stationary. To demonstrate these 
results, we introduce a powerful polynomial representation of commuting Majorana Hamiltonians. 
Remarkably, the physical properties of the topologically-ordered state are encoded in an algebraic 
variety, defined by the common zeros of a set of polynomials over a finite field. This provides a 
“geometric” framework for the emergence of topological order. 


Topological phases of matter are remarkable quantum 
states with quantized properties that are stable under 
local perturbations and can only be measured by non¬ 
local observables [2]. The most celebrated example is 
the fractional quantum Hall state, discovered more than 
thirty years ago [3]. The field of topological matter has 
now become an exciting research frontier at the cross¬ 
roads between theoretical physics, quantum information 
and material science. 

Our theoretical understanding of topological matter is 
largely built on topological quantum field theory (TQFT) 
[4]. In this framework, the action of quantum fields in 
a space-time manifold is independent of its metric, but 
depends crucially on its topology. Canonical quantiza¬ 
tion of these fields in a multiply-connected space yields 
a finite-dimensional Hilbert space, describing the degen¬ 
erate ground-states of topological matter. Wilson lines 
describe the world-lines of quasi-particle excitations, and 
the expectation value of “knotted” Wilson lines deter¬ 
mines the quasi-particle braiding statistics. A hallmark 
of topologically-ordered states in two dimensions is the 
presence of mobile quasi-particles with fractional statis¬ 
tics, or anyons [5]. 

Exactly solvable models often provide ideal play¬ 
grounds and valuable insights in theoretical studies of 
topological phases. In the past, a wide array of non- 
chiral topological phases in two dimensions have been 
obtained in spin models [6, 7], whose universal proper¬ 
ties are captured by topological quantum field theories. 
Recently, an exotic quantum phase with extensive topo¬ 
logical ground state degeneracy was discovered by Haah 
in three-dimensional (3D) spin models [8]. A remarkable 
property of this phase is that all topological excitations 
are strietly loealized in space, a feature which lies beyond 
the paradigm of topological quantum field theory. 

In this work, we introduce a wide range of 
translationally-invariant, solvable Hamiltonians of inter¬ 
acting Majorana fermions that exhibit a new kind topo¬ 
logical quantum order. These models have extensive 
topological degeneracy and a hierarchy of topological ex¬ 
citations that are only free to move within sub-manifolds 
of the full lattice. In one particular Hamiltonian in 


d = 3 spatial dimensions, the fundamental excitations 
are strictly localized, while composites of these excita¬ 
tions are free to move along one- and two-dimensional 
surfaces. The fundamental excitations are termed “frac- 
tons”, as they behave as fractions of a mobile particle. 

To systematically search for these models, compute 
their ground-state degeneracy on a d-dimensional torus 
and study their excitations, we introduce a purely alge¬ 
braic description of commuting Majorana Hamiltonians. 
We demonstrate that on a d-dimensional lattice with a 
two-site basis and a single interaction term per unit cell, 
an ideal Majorana Hamiltonian generally exhibits exten¬ 
sive topological degeneracy. We emphasize that each of 
our models may be written in terms of complex fermions 
by choosing appropriate pairings of Majorana fermions 
over the entire lattice. Our models also admit a loeal 
mapping to a boson model with identical topological de¬ 
generacy and a similar dimensional hierarchy of excita¬ 
tions, after projecting out half of the Hilbert space. We 
note that one of our models has similar phenomenology 
to a spin model studied in ref. [10, II]. 

Our approach to studying ideal Majorana Hamiltoni¬ 
ans provides a novel geometric framework for topological 
order, beyond topological quantum field theory. Remark¬ 
ably, a commuting Majorana Hamiltonian on a torus 
specifies an algebraie variety - defined as the common 
zeros of a collection of polynomials over a finite field - 
that encodes all physical properties of the topologically- 
ordered state. While a TQFT assigns a ground-state 
sector to an isotopy class of smooth, closed curves on a 
manifold, our models associate ground-state sectors with 
curves based on finer equivalence relations, resulting in 
extensive topological degeneracy in dimensions d > 3. 
We emphasize that our models are distinct from the ex¬ 
otic phase realized by Haah’s code [8], due to the presence 
of mobile composite topological excitations. 


I. OVERVIEW 

Due to the length of this paper, we begin with a de¬ 
tailed summary of our findings. We consider exactly 
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solvable Hamiltonians of interacting Majorana fermions 
that realize exotic forms of topological order. On a d- 
dimensional lattice with a basis, these Hamiltonians will 
be the sum of a single type of local operator over all 
lattice sites 


H = -J20m ( 1 ) 

m 

SO that all operators mutually commute and square to 
the identity, i.e.. 


[O^,On]=0, (2) 

(0„)2 = +1. (3) 

The operator On is required to be a product of an even 
number of Majorana fermions, so that the fermion parity 
of the entire system is conserved. A ground state |T) of 
(1) will satisfy the constraint that 

On^l^> = l^), (4) 


for all m. 

In Section H, we introduce a purely algebraic approach 
to systematically search for topological order in commut¬ 
ing Majorana Hamiltonians (1). A similar approach has 
been used previously to study topological order in com¬ 
muting Pauli Hamiltonians [9]. We represent the opera¬ 
tor O appearing in (1) as a set of Laurent polynomials 
over the field F 2 , which consists of two elements {0,1} 
with Z 2 addition and multiplication. We derive a math¬ 
ematical condition for a set of such polynomials to repre¬ 
sent a commuting Majorana Hamiltonian with topologi¬ 
cal order. This polynomial representation enables us to 
analytically determine the topological ground state de¬ 
generacy on a d-dimensional torus and deduce properties 
of topological excitations using algebraic methods. 

Using this polynomial approach, we demonstrate the 
following remarkable results. First, a topologically- 
ordered commuting Majorana Hamiltonian on a lattice 
with a two-site basis may be entirely specified by a sin¬ 
gle polynomial over F 2 . The ground state degeneracy for 
such a Hamiltonian on a d-dimensional torus of size L, 
which we denote by Dq, will generally take the asymp¬ 
totic form: 

log2L>0 = cL‘'-2 + 0(L'^-3) (5) 

for some constant c. We perform an exhaustive analy¬ 
sis and discover a class of commuting Majorana fermion 
models on a three-dimensional lattice with a two-site ba¬ 
sis, which exhibit extensive topological degeneracy of the 
form (5) with d = 3. 

Remarkably, despite being translationally invariant, 
our models admit fundamental point-like excitations that 
are strictly localized in space, and cannot move with¬ 
out paying a finite energy cost to create additional ex¬ 
citations. Composites of these fundamental excitations, 
however, are topological excitations that are free to move 
within sub-manifolds of the d-dimensional lattice. We 



FIG. 1. Majorana Cubic Model: The Majorana cubic 
model is defined on a cubic lattice, as in (a), with a single 
Majorana fermion per lattice site (colored red). The operator 
On is the product of the 8 Majorana fermions at the vertices of 
a cube. The Hamiltonian is a sum of these local operators over 
every other cube (colored blue) in a checkerboard pattern. As 
any pair of operators either share exactly one edge or none, all 
operators mutually commute. We choose to label the cubic 
operators A, H, C, and D as shown in (b). Acting with a 
single Majorana operator 7 ^ creates these four excitations. 


term these fundamental excitations that behave as frac¬ 
tions of mobile particles, “fractons.” Furthermore, we re¬ 
fer to bound states of fractons that can only move freely 
along an n-dimensional manifold as “dimension-n” par¬ 
ticles. In particular, a dimension-2 particle can be an 
anyon with well-defined fractional statistics. 

To motivate further study of ideal Majorana Hamilto¬ 
nians, we now describe in detail the phenomenology of 
fracton excitations and their composites in the simplest 
of our models, the Majorana cubic model. As shown in 
Figure 1(a), here the operator On is the product of the 
eight Majorana fermions at the vertices of a cube. The 
Hamiltonian is simply the sum of these operators over a 
face-centered-cubic (fee) array of cubes, forming a three- 
dimensional checkerboard. Since adjacent cubes share a 
common edge with two vertices, operators On on differ¬ 
ent cubes are mutually commuting, and their common 
eigenstate defines the ground state. For convenience in 
later analysis, we choose to identify four species of cube 
operators - A, H, C, and B - as shown in Figure 1(b). 

A fundamental excitation in the Majorana cubic model 
is obtained when the eigenvalue of a cube operator On 
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Excitation 

Type 

Statistics 

Operator 

ABCD 

Majorana 

Fermion 

7 

AA,BB, 
CC, DD 

Dim.-2 Anyon 

Boson 

Pair of Adjacent 
Wilson Lines 

AB, AC, 
AD, BC, 
BD, CD 

Dim.-l Particle 


Single Wilson 
Line 

A, B, C, D 

Fracton 


Membrane 


TABLE L Hierarchy of excitations in the Majorana cubic 
model. The fundamental cube excitation is a fracton, while 
two-fracton bound-states can behave as particles that are ei¬ 
ther free to move along one- or two-dimensional surfaces. The 
operator that creates each type of excitation is indicated. 

is flipped. The product of On over all cubes of a single 
type (A, C, or D) is equal to the fermion parity T of 
the entire system and is fixed. 

T=l[0,= ll0p=l[0p=ll0p, (6) 

p^A p£B pEC p£D 

Therefore, a single cube-flip excitation cannot be created 
alone, and is a topological excitation. Remarkably, the 
fundamental cube excitation in this model is completely 
immobile^ as we observe through the following physical 
argument. In the cubic model, acting on the ground-state 
with a single Majorana fermion flips the eigenvalues of 
four adjacent cube operators, as shown in Figure 1(b). 
This four-cube excitation may trivially move by acting 
with a Majorana bilinear. If the fundamental cube ex¬ 
citation were mobile, then it would be possible to move 
it in any arbitrary direction, as the cube operator it¬ 
self preserves all lattice symmetries. In this case, the 
cube excitation would have well-defined (fermion or bo¬ 
son) statistics, and a four-cube bound-state could never 
be a fermion. Therefore, it must be the case that the 
fundamental cube excitation is frozen. A rigorous proof 
of the immobility of the fundamental excitation is given 
in Section IV using the polynomial representation of the 
ideal Majorana Hamiltonian. 

We now analyze the fracton bound-states in the Majo¬ 
rana cubic model in detail, along with the mutual statis¬ 
tics of the excitations. Using the labeling of the cube 
operators shown in Figure 1(b), we find the hierarchy of 
quasiparticles shown in Table I in the Majorana cubic 
model. The fundamental fracton excitation appears at 
the corners of membrane-like operators and may only be 
created in groups of four. Two-fracton bound-states can 
form dimension-1 particles or dimension-2 anyons. Re¬ 




FIG. 2. Dimension-1 Particle: Excitations (colored) may 
be created by acting with Wilson line operators. In (a), a 
straight Wilson line creates pairs of dimension-1 particles at 
the endpoints. The dimension-1 particle may hop freely in the 
direction of the Wilson line, by acting with Majorana bilinear 
terms. Remarkably, the dimension-1 particle cannot hop in 
any other direction without creating additional excitations. 
Introducing a “corner” in the Wilson line, as in (b), creates 
an additional topological excitation localized at the corner. 

markably, a dimension-2 anyon has tt mutual statistics 
with a fracton lying in its plane of motion. As a result, 
while the fracton is immobile, its presence may be de¬ 
tected by a braiding experiment. Furthermore, the exact 
location of a single fracton within a finite volume V may 
be determined by braiding dimension-2 anyons in the 
three mutually orthogonal planes around the boundary 
dV. In this way, the exact quasiparticle content within 
V is effectively encoded “holographically”, and may be 
determined by ^ 0(i) braiding experiments, where i is 
the linear size of a box bounding V. 

We now proceed to explore the hierarchy of excitations 
in detail. 

Dimension-1 Particle: The dimension-1 particle 
may be created by acting with a single Wilson line oper¬ 
ator, defined by the product of the Majorana operators 
along a straight path i. Up to an overall pre-factor of 
±1, ±i, we write the Wilson line operator as 

We^Yljn. (7) 

nG£ 

As shown in Figure 2(a), the straight Wilson line anti¬ 
commutes with two cube operators at each of its end¬ 
points; the two cube excitations at a given endpoint are 
of different types. As a result, Wi creates pairs of exci- 
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tat ions of the form AB^ AC ^ AD^ BC, BD^ or CD. Re¬ 
markably, these two-fracton bound-states are only free to 
move along a line, by simply extending the Wilson line 
operator Wi by acting with a Majorana bilinear along the 
path i. If we try to move this two-fracton bound-state in 
a plane, we must introduce a corner in the Wilson line, 
which localizes an additional topological excitation at the 
corner, as shown in Figure 2(b); the excitation cannot be 
removed by the action of any local operator. As the pat¬ 
tern of excitations produced by a Wilson line Wi is sen¬ 
sitive to the geometry of i, the two-fracton bound-states 
AB, AC^ AD^ BC^ BD^ and CD are restricted to move 
along a line and behave as dimension-1 particles. We em¬ 
phasize that they cannot move in a higher-dimensional 
space without creating additional cube excitations. 

Dimension-2 Anyon: Acting with a pair of adjacent 
Wilson lines and along parallel paths I and 

respectively, also creates a pair of two-fracton bound- 
states localized at the ends, as shown in Figure 3(a). At 
each end of the path, however, the operator 
now creates pairs of cube excitations of the same type 
(AA, BB^ CC or DD). These two-fracton bound-states, 
where each fracton is of the same type, are allowed to 
move freely in the two-dimensional plane orthogonal to 
the shortest line segment connecting the two paths i and 
C without creating additional excitations; this is shown 
in Figure 3(b). We note that detailed geometric fea¬ 
tures of a single Wilson line, such as the presence of 
sharp corners, determine the pattern of excitations cre¬ 
ated from the ground-state. However, when acting with 
an appropriate pair of adjacent Wilson lines, the excita¬ 
tions created at the sharp corners may be annihilated. 
Therefore, a pair of adjacent Wilson lines may be de¬ 
formed in the plane with no energy cost. We conclude 
that the AA, BB^ CC and DD two-fracton bound-states 
are dimension-2 anyons. 

Braiding a dimension-2 anyon around a closed loop in 
the plane is equivalent to acting with the product of cube 
operators within the two-dimensional region enclosed by 
the loop; this is shown for a particular choice of loop in 
Figure 3(c). As a result, braiding a dimension-2 anyon 
around a closed loop enclosing a single fracton in the 
plane produces an overall minus sign. The ability to de¬ 
tect a fracton with a dimension-2 anyon produces non¬ 
trivial mutual statistics between the dimension-2 anyon 
and other particles in the excitation spectrum of the Ma¬ 
jorana cubic model. First, a dimension-2 anyon has tt 
mutual statistics with any dimension-1 particle in the 
same plane, as braiding the dimension-2 anyon in a closed 
loop will only detect one of the two fractons that make up 
the dimension-1 particle. Furthermore, the dimension-2 
anyon has tt mutual statistics with dimension-2 anyons 
that are free to move in adjacent, parallel planes. 

Fractons and Membrane Operators: Acting with 
Majorana operators on a flat, two-dimensional membrane 
S creates fracton excitations at the corners of the bound¬ 
ary of S, as shown in Figure 4. We write the membrane 



(c) 


FIG. 3. Dimension-2 Anyon: Acting with two adjacent 
Wilson line operators Wi and W 2 creates pairs of excita¬ 
tions at the endpoints of the same type (AA, BB, CC or 
DD). These two-fracton excitations are free to move in a 
two-dimensional plane orthogonal to the shortest line seg¬ 
ment connecting the pair of Wilson lines. Furthermore, in 
(b) we may detect a fracton (colored blue) by braiding a 
dimension-2 anyon around a closed loop enclosing the frac¬ 
ton. As the braiding operator, a pair of closed Wilson line 
operators W 1 W 2 , is equal to the product of the enclosed cube 
operators as shown above. Therefore, the braiding produces 
an overall minus sign if an odd number of fractons are en¬ 
closed. 


operator up to an overall pre-factor of ±1, as 

(8) 

For a rectangular membrane in the x-y plane, the bound¬ 
ary 511 is a closed, rectangular loop with dimensions 
and ly. We note that if and ly are both even, then 
the fracton excitations created at the corners of will 
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FIG. 4. Membrane Operator &; Practon Excitations: 

Acting with a product of Majorana operators on a surface E 
creates localized excitations at the corners of the boundary 
^E as shown above. 


all be of the same type. Alternatively, if is odd and ^y 
is even, then the pairs of fracton excitations separated in 
the ^-direction will be of the same type, while fractons 
separated in the x-direction will be distinct. 

Extensive Topological Degeneracy: Using the al¬ 
gebraic representation of the Majorana cubic model, we 
compute its ground-state degeneracy Dq to be 

log2 Do = 31/ — 3 (9) 

on an 1/ X 1/ X 1/ three-torus, with periodic boundary 
conditions imposed in the x, and ^ directions, with 
each cube having unit side length. Pairs of string-like 
Wilson loop operators wrapping non-trivial cycles of the 
torus - corresponding to tunneling dimension-2 anyons - 
distinguish the ground-state sectors. As the number of 
distinct dimension-2 anyons grows linearly with system 
size, the ground-state degeneracy is necessarily extensive. 

We emphasize that the algebraic approach allows us 
to systematically search for topologically-ordered, ideal 
Majorana Hamiltonians, rigorously characterize the na¬ 
ture of excitations, and calculate the ground-state de¬ 
generacy in a wide range of Majorana models using tech¬ 
niques in algebraic geometry. As a result, the next two 
sections of this work introduce and focus on the polyno¬ 
mial representation of ideal Majorana Hamiltonians and 
draw broad conclusions based on this representation. In 
Section HI, we present the 6 distinct three-dimensional 
Majorana models with nearest-neighbor interactions that 
are topologically-ordered. In particular, one of our mod¬ 
els, which may naturally be written in terms of complex 
fermions on an fee lattice, has a fundamental excitation 
that may only freely move along a line in the (1,1,1) di¬ 
rection. 

We conclude, in Section IV, with a proof of the pres¬ 
ence of fractons in the Majorana cubic model, and briefly 
outline the phenomenology of excitations in the remain¬ 
ing models. 


II. TOPOLOGICAL ORDER IN COMMUTING 
MAJORANA HAMILTONIANS 

In this section, we introduce a representation of the op¬ 
erators in the ideal Majorana Hamiltonian (1) as a vector 
of Laurent polynomials over the finite field F 2 . The alge¬ 
braic representation provides an important starting point 
for studying and classifying Majorana Hamiltonians. We 
demonstrate that the following conditions, that 

(i) All operators in the ideal Hamiltonian mutually 
commute, and 

(ii) Degenerate ground-states of the Hamiltonian are lo¬ 
cally indistinguishable 

may be phrased entirely in the polynomial representa¬ 
tion. The ground-state degeneracy of an ideal Majorana 
Hamiltonian (1) on the torus can be computed as the 
dimension of a quotient ring [9] . 

We demonstrate that an ideal Majorana Hamiltonian 
obeying (i) and (ii) on a lattice with a two-site basis 
and a single interaction term per unit cell may be spec¬ 
ified by a single polynomial over F 2 . We use this result 
to systematically search for and characterize commuting 
Majorana Hamiltonians. In d = 3 dimensions, we find 6 
distinct, non-trivial models with nearest-neighbor inter¬ 
actions, extensive topological degeneracy, and a dimen¬ 
sional hierarchy of excitations. 

A. Algebraic Representation 

To study commuting Majorana Hamiltonians, we rep¬ 
resent the operator O appearing in Eq. (1) as a polyno¬ 
mial over the field F 2 . A similar mapping has been in¬ 
troduced in the context of Pauli Hamiltonians [9] . Con¬ 
sider a d-dimensional lattice with translation operators 
{ti,... ,td} and an n-site unit cell. We restrict n to be 
an even integer so that there is a well-defined number 
of complex fermions per lattice site. We label the Ma¬ 
jorana fermions within the unit cell at the origin as jj 
for j = 1, 2,..., n. All other Majorana fermions on the 
lattice are obtained by acting with translation operators. 

Any Hermitian operator acting on this lattice may be 
written as a sum of products of Majorana operators. For¬ 
mally, we may write a summand O as 

n 

n ■•■c(10) 

with Ui e Z and Cj(ni,... G {0,1}. For simplicity, 
we have omitted the prefactor in the expression 

for O, which plays no role in our analysis. We introduce a 
purely algebraic representation of this operator by noting 
that any product of translation operators may be writ¬ 
ten as a monomial, e.g. In 

this way, the action of the translation group is naturally 
represented by monomial multiplication. 
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Recall that distinct Majorana fermions anti-commute 
and that each Majorana operator squares to the identity. 
Therefore at each site within a unit cell, the identity 1 
and 7 under multiplication form the group Z 2 , with the 
two operators represented by the group elements 0 and 1 , 
respectively. In this representation, the operator equal¬ 
ity 7 ^ = 1 maps to the Z 2 group addition 1 + 1 = 0 . 
This simple algebra of Majorana fermions allows us to 
write any product of Majorana operators as the sum of 
monomials - representing the location of each Majorana 
operator via the action of the translation group - with 
Z 2 coefficients. As an example, consider a lattice with 
a single site per unit cell, and the Majorana operator 7 
at the origin. A Majorana bilinear admits the following 
polynomial representation: 


^ l + ( 11 ) 

In this notation, operator multiplication corresponds to 
polynomial addition with Z 2 coefficients. 

For the general case of a unit cell with n sites,we rep¬ 
resent a product of Majorana operators as a vector of 
polynomials over F 2 , with the j-th entry of the vector 
representing the action of the translation group on 7 ^, 
the j-th Majorana fermions in the unit cell at the origin. 
For example, the operator (10) may be written as 




{rii} 


y,nd 


ci(ni 

C2(ni 


rid) \ 

rid) 


( 12 ) 


\C2n(ni, ...,nd)J 


Adopting the terminology in Ref. [9] , we refer to S as the 
“stabilizer map” for the remainder of this work. 

To illustrate the algebraic representation of operators 
in commuting Majorana Hamiltonians, we present a con¬ 
crete example. Consider the Majorana plaquette model 
in Ref. [1], which is defined on a two-dimensional honey¬ 
comb lattice with one Majorana fermion per site and a 
Hamiltonian of the form (1) where Op is the product of 
the six Majorana fermions at the vertex of a hexagonal 
plaquette p. We show a single hexagonal plaquette on the 
lattice in Figure 5(a), along with the Majorana fermions 
7 a and 75 within the two-site unit cell. The correspond¬ 
ing stabilizer map S{x,y) for the six-Majorana operator 
is given by: 



Consider a honeycomb lattice with a single Majorana fermion 
on each lattice site. We define an operator Op as the product 
of the six Majorana fermions on the vertices of a hexagonal 
plaquette p, as shown in (a). The colored plaquettes in (b) 
correspond to the three distinct bosonic excitations (A, H, or 
C) that may each be created in pairs by acting with Wilson 
line operators. 


creates an excitation. We use a polynomial to record 
the locations of all excitations in the state IF |T); each 
location is labeled by the translation vector connecting 
it to the origin. Specifically, for a Hamiltonian with 
stabilizer map 5'(xi,--- ^Xd) and an arbitrary operator 
IF with a polynomial representation P(W) of the form 
( 12 ), we define the “excitation map” E{xi^... ^Xd) so 
that E{xi ,..., Xd) • P{W) G F 2 [xf^, • • • ,x^^] describes 
the excitations created by IF. In the Supplemental Ma¬ 
terial [12], we demonstrate that E is simply given from 
the stabilizer map as follows: 


E{xi,...,Xd) = ,Xd) (14) 


where S{xi, ...,Xd) = [^(a;!, .. .,Xd)]'^. 


As an example, the excitation map for the Majorana 
plaquette model is given by E{x^ y) = (1 + x + p, 1 + x + 
xy ). Below, we show the action of the operator 7 ^ at the 
origin in the Majorana plaquette model, which creates 
three adjacent excitations as specified by the red points. 
The locations of the excitations are obtained by perform¬ 
ing the matrix multiplication of E with the polynomial 


representation 



of 7 a: 


/l+x+p\ 

S{x,y) = \ . (13) 

ylPxPxyJ 

Here, we adopt the notation that y = y~^^ x = x~^. 
As shown in Ref. [1], this Hamiltonian exhibits a novel 
form of Z 2 topological order with fermion parity-graded 
excitations and exact anyon permutation symmetries. 

Next, we consider the action of an arbitrary operator 
IF on the ground state |F) of the commuting Majorana 
Hamiltonian. When IF anticommutes with an operator 
On in the Hamiltonian, it flips its eigenvalue and thus 


E{x,y)- =l + x + y. (15) 

Therefore, the action of 7 ^ may be represented by the 
polynomial IPxPy^ labeling the locations of the flipped 
plaquettes; here, the plaquette operator corresponding to 
the origin (i.e. the location “ 1 ”) is to the right of 7 ^, as 
can be seen from its polynomial representation (13). 

A dictionary that summarizes the relationship between 
Majorana operators and polynomials is given in Table H. 
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Operator 

Polynomial 

Majorana Fermion 

Vector over F 2 



[j = 1 ,..., n for each 
site in the unit cell] 

[n-dimensional unit vector 
with entry equal to 1] 

Translation 

Monomial Multiplication 



Multiplication 

Addition in F 2 [xi, • • • , Xd] 

ij ■ tilt 

Cj + {xkY eg 


TABLE IL Summary of the polynomial representation of Ma- 
jorana operators. An arbitrary operator in d spatial dimen¬ 
sions, written as the product of Majorana fermions, may be 
represented as a vector with entries in the (Laurent) polyno¬ 
mial ring F 2 [xf ^, • • • , X . 


B. Topological Order and Ground-State 
Degeneracy in the Algebraic Representation 

The polynomial representation of Majorana operators 
serves as a starting point for constructing commuting 
Majorana Hamiltonians that exhibit topological orders. 
As we demonstrate in the Supplemental Material [12], for 
a translationally invariant Majorana Hamiltonian with a 
single operator per lattice site, all operators mutually 
commute if and only if its stabilizer map 5'(xi, • • • ,Xd) 
satisfies the condition 


S{xi,...,Xd)- S{xi,...,Xd) =0. (16) 

More generally, if the Hamiltonian contains multiple op¬ 
erators per lattice site then we may define a set 

of stabilizer maps for each type of operator {5'i}, so that 
the condition ..., Xd) • ..., Xd) = 0 for all i, 

j, guarantees that all terms in the Hamiltonian commute. 

We next formulate a necessary and sufficient alge¬ 
braic condition for topological order in commuting Ma¬ 
jorana Hamiltonians, which requires that any degener¬ 
ate ground-states of a topologically-ordered Hamiltonian 
cannot be distinguished by local operators. The local in- 
distinguishability is equivalent to the condition that, for 
any local operator Mi 

ncsM^ncs = c(M^)nGs (i7) 

where Hgs is the projector onto a ground-state sector and 
c{Mi) is a constant that only depends on the operator. 
For our case, consider an operator M/ that is the product 


of Majorana operators, and P(M^), the polynomial rep¬ 
resentation of Mi. If Mi anti-commutes with any term in 
the Hamiltonian, then Mi creates excitations when acting 
on the ground-state, and we have Hgs Mi Hgs = 0 - If Mi 
commutes with the Hamiltonian, then P{Mi) G kerP, 
as Mi creates no excitations. In this case, the condi¬ 
tion Hgs Mi Hgs = c(M^)nGs is guaranteed if Mi may 
be written as a product of operators already appearing 
in the Hamiltonian. More generally, any local operator 
M that commutes with the Hamiltonian then takes the 
form: 

M = ^Mi (18) 

i 

where each term Mi is the product of operators already 
appearing in the Hamiltonian. This condition is neces¬ 
sary for distinct ground-state sectors to be locally indis¬ 
tinguishable. 

In our polynomial representation, we enforce the con¬ 
dition (17) by requiring that the stabilizer and excitation 
maps satisfy the following condition on an infinite lattice 

keTE = imS. (19) 

Recall that the image of S is the set of all polynomial 
linear combinations of 5'(xi,..., x^^), taking the form of 

^x^^ ■■■x2‘'S{xi,...,Xd), (20) 

{ni} 

and representing all operators that can be written as a 
product of the commuting operators appearing in the 
Hamiltonian. On the other hand, the kernel of the ex¬ 
citation map E is the set of all operators that do not 
create any excitation when acting on the ground state. 
The above algebraic condition (19) for topological order 
is thus equivalent to the statement that any operator 
that creates no excitations on a ground state on an in¬ 
finite lattice is necessarily a product of operators {On} 
already appearing the Hamiltonian. In other words, there 
are no non-trivial, locally conserved quantities, and any 
degenerate ground-states of the Hamiltonian are locally 
indistinguishable. In summary, imposing the commuta¬ 
tivity (16) and local indistinguishability (19) conditions 
on a stabilizer map produces an ideal Majorana Hamil¬ 
tonian with topological order. 

We may compute the ground-state degeneracy of an 
ideal Majorana Hamiltonian in the polynomial represen¬ 
tation via constraint-counting. A lattice with 2M Ma¬ 
jorana fermions defines a 2^-dimensional Hilbert space. 
On the torus, however, fixing the eigenvalues of the com¬ 
muting operators in the ideal Majorana Hamiltonian only 
imposes M — k multiplicatively independent constraints, 
since the product of certain operators appearing in the 
Hamiltonian will yield the identity. The ground-state de¬ 
generacy is simply given by the space of states satisfying 
the constraints, which is precisely 2^/2^“^ = 2^. As 
each ideal Majorana Hamiltonian in this work consists of 
exactly one term for each pair of Majorana modes, we see 
















that k is directly equal to the number of constraints on 
the commuting operators appearing in the Hamiltonian. 

For example, in the Major ana plaquette model, we may 
group the plaquette operators {Op} into three types (H, 
H, and C) as shown in Figure 5(b). On the torus, the 
product of the H, H, and C-type operators is identical 
and equal to the total fermion parity [1]. This yields the 
following two independent constraints: 

n Op n =n n (^i) 

peA peB peB pec 

and produces a 2 ^-fold degenerate ground-state on the 
torus. These constraints may be compactly represented 
using polynomials labeling the locations of the H, B and 
C-type plaquettes. For example, the collection of all A 
plaquettes is captured by the polynomial 


PA = {l + xy + 



It is straightforward to expand pa to verify that the ex¬ 
ponents of the non-zero terms describe the positions of A 
plaquettes. Here, L specifies the periodic boundary con¬ 
ditions in the x and y directions, so that = 1 , = 1. 

Similarly, the collections of all plaquettes in B and C are 
encoded in ypA and xypA^ respectively. The constraints 
( 21 ) arise from the fact that {pA + yPA)S = 0 , using 
Eq. (13) and the boundary conditions. 

In terms of the stabilizer map, any multiplicative con¬ 
straint on the operators in the ideal Majorana Hamilto¬ 
nian on the torus is in one-to-one correspondence with a 
solution p of the equation p • S' = 0 , so that the polyno¬ 
mial p is an element of the kernel of S. Therefore, the 
number of independent relations is given by 


k = dimF 2 [ker(S)]. (23) 

We rewrite the expression (23) in a more convenient 
form for calculations that will also allow us to make gen¬ 
eral statements about the scaling behavior of the ground- 
state degeneracy with system size for an ideal Majorana 
Hamiltonian of the form ( 1 ). As proven in Corollary 4.5 
in Ref. [9], Eq. (23) is equivalent to the dimension of the 
following quotient ring: 


k = log 2 D = dimF 2 


/F2[xi,-- - ,Xd\\ 

V HS) + bL )' 


(24) 


Here, I{S) the ideal generated by the stabilizer map; if 
= {si,, S 2 n) then I{S) is the space of polynomials 
in F 2 [a:y,... ,xp\ obtained as a linear combination of 




I{S) 



Ci G F2Ff\ . 



(25) 


We will denote the ideal generated by a set {si,..., 
by (si,..., Sn)‘ Eurthermore, we define the ideal bi, = 


(xf — 1 ,..., — 1 ). As the quotient space identifies the 

zero element in F 2 [xf^, ... with the generators of 

I{S) + we observe that the ideal bi, is used to enforce 
the periodic boundary conditions on a d-dimensional 
torus with side-length L. 

We emphasize that the ideal I{S) is the space of 
excitations that can be created through the action of any 
operator on the ground-state. Therefore, the expression 
(24) may be physically interpreted as counting certain 
superselection sectors of the ideal Majorana Hamilto¬ 
nian. Any p G ¥ 2 [xf^, • • • ,corresponds to a virtual 
eigenstate of the Hamiltonian with excitations at the 
locations specified by the polynomial p. Certain states, 
however, cannot be created by acting with an operator 
on a ground-state |T) due to the k constraints on the 
commuting operators. Eor example, in the Majorana 
plaquette model, it is impossible to obtain a state with 
a single plaquette excitation by acting on the ground- 
state, since the products of A, B and C plaquettes 
must satisfy (21). As I{S)/{I{S) fl bi,) is the set of 
excitations that can be created by the action of operators 
on the ground-state for a finite system, the quotient 
space (F 2 [a;f\--- ,a;^V^’L)/U('S')/(-f(S')nbL)) = 

F 2 [xi,-'- ,Xd]/{bL¥l{S)) is the set of virtual eigen¬ 
states of the Hamiltonian that cannot be deformed into 
each other through the action of any local operator. For 
the Majorana plaquette model, this quotient space is 

_ ^2[x,y] _^ ^2 (26) 

{1 + x + y, X + y + xy, x^ - 1, y^ - 1) ^ 

when L mod 3 = 0 so that there are an equal number of 
A, H, and C plaquettes. In this case, the trivial vacuum 
(0) and a state with a single plaquette excitation (1) on 
A, H, or C correspond to the four superselection sectors 
in the quotient ring. 

The expression for the ground-state degeneracy (24) 
is convenient as the dimension of a quotient ring may 
be computed using algebraic techniques. Most often, we 
will determine a Grobner basis for the ideal I{S) + bi, 
in order to determine membership in the quotient ring. 
For a polynomial ring R, we may define a total mono¬ 
mial ordering (e.g. lexicographic order with xi y X 2 ^ 
... ^ Xd); we denote the leading monomial in a poly¬ 
nomial h e R AS LM(h) with respect to this ordering. 
Given an ideal I = (si,--- ,s^) of a polynomial ring, 
there exists a canonical choice of generators for the ideal, 
known as the Grobner basis { 5 ^ 1 , •• • , with the prop¬ 
erty that for any / G /, LM(/) G (LM(^i),... ,LM(^n)), 
i.e. any element of the ideal has a leading term con¬ 
tained in the ideal generated by the leading terms of the 
Grobner basis. As a result, the dimension of the quotient 
ring dim[R//] is merely given by the number of monomi¬ 
als that are smaller (in the monomial ordering) than all 
of the leading terms in the Grobner basis. This is because 
any polynomial p C R may be reduced by the Grobner 
basis until the leading term of the reduced polynomial 
satisfies LM(pred) < LM(^^) for all i = 1,... ,n. There¬ 
fore, each monomial m satisfying m < LM(^^) for all i 
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Algebraic Expression 

Physical Interpretation 

S(xi, ..,Xd) • S(xi,..,xd) = 0 

Commutativity condition, 
that all operators {On} 
appearing in the 

Hamiltonian mutually 
commute. 

im(S') 

Set of operators that may 
be written as the product of 
commuting operators {On} 
in the Hamiltonian. 

ker(E) 

Set of operators that create 
no excitations when acting 
on the ground-state T). 

k — dimF2 [ker(AS')] 

The number of independent 
relations among the 
commuting operators in the 
Hamiltonian, when placed 
on a torus. The ground 
state degeneracy D = 2^. 

ve¥2[xf,...,x^] 

A configuration of 
excitations, specified by 
the locations of operators 
{On} with flipped 
eigenvalue —1. 

q € I{S) 

A configuration of 
excitations that may be 
created by acting with an 
operator on the 
ground-state T). 


TABLE III. Dictionary of various algebraic quantities and 
their physical interpretation in the context of a commuting 
Majorana Hamiltonian. 


corresponds to a unique representative of the quotient 
ring Rjl. 

We note that calculations of the ground-state degen¬ 
eracy for any commuting Majorana Hamiltonians pre¬ 
sented in this work are done by determining a Grobner 
basis for the ideal 1[S) + bi,. In this way, the calcula¬ 
tion of the degeneracy reduces to counting points in an 
algebraic set. 

C. Unitary and Stable Equivalence 

The polynomial representation of the ideal Majorana 
Hamiltonian contains built-in redundancies, since we 
may re-define the unit cell or translation operators on 


the d-dimensional lattice. For the stabilizer map, the 
translation corresponds to multiplication of any entry of 
5'(xi,..., by a monomial. In this way, a stabilizer 
map 5'(xi,..., is only defined up to monomial mul¬ 
tiplication on each of its entries. Furthermore, for an 
ideal Majorana Hamiltonian with longer-range interac¬ 
tions, we may always enlarge the unit cell. As our focus 
will be on Majorana models with nearest-neighbor in¬ 
teractions, we neglect this redundancy in the stabilizer 
map. 

Equivalence relations, given by local unitary trans¬ 
formations on ideal Majorana Hamiltonians, may also 
be considered in the polynomial language. For in¬ 
stance, two ideal Majorana Hamiltonians, defined by sta¬ 
bilizer maps S{xi ,..., Xd) and 5"(xi,..., x^) are uni- 
tarily equivalent if there exists a matrix U such that 
5"(xi, ...,Xd) = U ' S{xi, ...,Xd) where U G 0 (n;F 2 ), 
an orthogonal ma trix over F 2 sa tisfying U^U = 1. This 
guarantees that if 5'(xi,..., Xd) • S{xi ^..., Xd) = 0 , then 
S'{xi ^..., Xd) • S'{xi ^..., Xd) = 0 as well. Finally, we 
take two stabilizer maps to be stably equivalent if we can 
obtain one from the other by attaching a trivial (dimer¬ 
ized) set of Majorana fermions. This is expressed as 
S{xi,... ,Xd)^ ^ 5'(xi,...,Xrf)^0 ( 0 , ••• , 0 , 1 , 1 ). 

III. EXTENSIVE TOPOLOGICAL 
DEGENERACY IN d > 3 

Using the commutativity (16) and local indistinguisha- 
bility (19) conditions, and the built-in redundancy in the 
polynomial description, we demonstrate in the Supple¬ 
mental Material [ 12 ], that an ideal Majorana Hamilto¬ 
nian defined on a d-dimensional lattice with a two-site 
basis is topologically-ordered if the stabilizer map may 
be written in the following form, after multiplying each 
entry by appropriate monomials: 

,Xd)\ 

S = _ (27) 

\f{xi,--- ,Xd)J 

where f{xi, ■ ■ ■ ,Xd) G and / and / are 

co-prime, i.e., / and / have no common polynomial fac¬ 
tors. As a result, a topologically-ordered, ideal Majo¬ 
rana Hamiltonian with a two-site basis may be speci¬ 
fied by a single polynomial. For example, the stabilizer 
map for the Maj orana p laquette model takes the form 
> 5 ^ = {f{x, y),x- f{x, y) ) with f{x, y) = l+ x + y. 

The dimension of the quotient ring (24) scales as the 
dimension of the space of the zeros of the ideal I{S) over 
the field extension ¥ 2 m when L = 2 "^ — 1. As a result, for 
an ideal Majorana Hamiltonian ( 1 ) with a two-site unit 
cell, the space of solutions to 

f{xi,...,Xd) = 0 , f{xi,...,Xd) =0 (28) 

generally defines an (d — 2)-dimensional variety, so that 
the ground-state degeneracy scales on the d-dimensional 
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h{x,y,z) 

f 2 {x,y,z) 

h{x,y,z) 

f 4 ,{x,y,z) 

Mx,y,z) 

h{x,y,z) 

/ 

1 + X + y + z 

1 + z-^xy 
-hyz + xz 

1 + X + y 
-Ayz + xz 

l-\-y A- z 
-\-xy + yz xz 

1 X y z 
-\-xy + yz xz 

1 + x + y 
+z -A yz 


■1 ' ~ ~A — 2 

A 7 - - t'' ! 

r I [ '1 

- ~ ^ 

r ' ' 1 

^Ax ! 



y% ! 

loga D 

3L-2 

6L-6 {L = 3n) 

0 {L ^ 3n) 

4L/3 (L = 6 • 2’") 
8L/5 h = 5 • 2”) 

4L - 4 (L = 2n) 

2L - 1 (L = 2n+ 1) 

2L-2 (L = 2^”+^ - 1) 
2L - 4 (L = 2^'* - 1) 

4L/3 (L = 6 • 2’") 
8L/5 h = 5 • 2”) 


TABLE IV. We find 7 distinct, topologically-ordered ideal Majorana Hamiltonians with nearest-neighbor interactions on a 
lattice with a two-site unit cell in d = 3 spatial dimensions. The first model fo{x,y,z) = 1 + y + z (not shown) is a trivial 
stack of two-dimensional Majorana plaquette models, considered in Ref. [1]. For the remaining 6 models, the action of the 
elementary operator O appearing in the ideal Majorana Hamiltonian is shown above as the product of the Majorana fermions 
on the indicated red dots. In the depiction of the Majorana cubic model fi{x, y, z), we have also shown the choice of translation 
vectors {tx,ty,tz} on the lattice, originating from one of the sites within the unit cell; to compute the ground-state degeneracy 
on an L X L X L torus, we impose periodic boundary conditions by requiring that = 1. The topological ground- 

state degeneracy (D) of each of these models is extensive. For models f 3 (x,y,z), / 5 (x,y,z), and f6{x,y,z), the ground-state 
degeneracy on the three-torus is a highly sensitive function of system size, and only the maximum value of the degeneracy is 
shown for the indicated choices of L. 


torus with side-length L as log 2 l^ = cL^~‘^ H-for some 

constant c. We emphasize that this produces a class of 
ideal Majorana models with extensive topological degen¬ 
eracy in d = 3 dimensions. Remarkably, while our models 
have a two-dimensional Hilbert space and a single inter¬ 
action term per lattice site, this only constrains the full 
Hilbert space up to extensive topological degeneracy. 

We have exhaustively searched for distinct, ideal Ma¬ 
jorana Hamiltonians with a two-site basis and nearest- 
neighbor interactions in d = 2 and d = 3 spatial dimen¬ 
sions. This is straightforward as the orthogonal group 
0 ( 2 ; F 2 ) = { 12 x 25 so that the space of local unitary 
transformations between these ideal Majorana Hamilto¬ 
nians is trivial. In d = 2 spatial dimensions, the only 
such Hamiltonian is the Majorana plaquette model with 

f{x,y) = l + x + y. (29) 

In d = 3 dimensions, however, we find 7 distinct Ma¬ 
jorana models with a two-site basis and nearest-neighbor 
interactions. The first model has the polynomial repre¬ 
sentation fo{x,y,z) = 1 +^ + 2 ) and is a trivial stack of 
two-dimensional Majorana plaquette models. The poly¬ 
nomial representations of the remaining models, along 
with their ground-state degeneracies on a torus of side- 
length L are shown in Table IV. For simplicity, we have 
imposed periodic boundary conditions by requiring that 
= ty = = 1 for the translation vectors {tx^ty^tz} 

shown in the representation of the Majorana cubic model 
fi{x,y,z) = 1 + X + ^ + 2 ) in Table IV. Each of the mod¬ 
els shown exhibits extensive topological degeneracy and 
admits at least one topological excitation that is free to 


move in a sub-manifold of the full lattice. 


IV. FRACTON EXCITATIONS AND 
DIMENSION-n ANYONS 

A remarkable feature of these Majorana models is the 
presence of fundamental excitations that are either per¬ 
fectly immobile or only free to move in a sub-manifold 
of the lattice; attempting to move these excitations by 
acting with any local operator will necessarily create ad¬ 
ditional excitations. A bound-state of these immobile ex¬ 
citations, however, forms a particle that can freely move 
along a higher-dimensional sub-manifold. 

The existence of a fracton fundamental excitation may 
be shown rigorously in the polynomial representation of 
the Majorana models. An element p G I{S) of the ideal 
defined by the stabilizer map corresponds to a set of ex¬ 
citations that may be created by acting on the ground- 
state. The fundamental excitation is mobile if and only if 
it is possible to create an isolated pair of such excitations. 
Therefore, an ideal Majorana model admits fracton exci¬ 
tations if the stabilizer ideal contains no binomial terms, 
i.e. 

l + x^^x^^ ---xy ^I{S) (30) 

for any G Z. 

We now apply the polynomial criterion for fracton ex¬ 
citations to the Majorana cubic model and to the model 
f^{x^y,z) = l-\-x^y-\-z-\-xy^yz-\- xz, both shown in 
Table IV. 
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A. Practons in the Majorana Cubic Model 

We consider the Majorana cubic model, specified by 
the single polynomial /(x, z) = so that the 

stabilizer map is given hy S = (^f{x^y^z)^ f{x^y^z)^ . 

We wish to prove that the ideal generated by the stabi¬ 
lizer map I{S) contains no binomial terms, so that the 
fundamental cube excitation is a fracton. This may be 
shown by considering the zero-locus of the ideal, i.e., the 
solutions to the zeros of the generators of the ideal: 

l + x^y + z = 0 (31) 

xyz xy ^ yz ^ xz = 0. (32) 

A polynomial p belong to I{S) only if p vanishes on the 
zero-locus of the ideal. Note that solutions to (31) take 
the form (x, y, z) = ( 1 , a, a), (a, 1 , a) or (a, a, 1 ), where 
a is an arbitrary element in the extension of F 2 . However, 
we see that the binomial 1 -h x^y^z^ vanishes on this 
space of solutions only ifn = m = ^ = 0 , in which case 
the binomial is zero. Therefore, we conclude that 

l + x^y^z^ ^ I{S). (33) 

As a result, there is no way to create the fundamental 
cube excitation in the Majorana cubic model in pairs. 
Therefore, the cube excitation is an immobile fracton; a 
single cube excitation cannot be moved without creating 
additional excitations. 


B. Dimension-1 Fundamental Excitations in 

Mx,y,z) 

Now, we consider the isotropic model f 5 {x, y,z) = 1 + 
x + y + z + xy + yz + xz, with stabilizer map defined 


by S{x,y,z) = [f^{x,y,z), xyz ■ f^{x,y,z)^. Prom the 

excitation map E{x^y^z) = S{x,y,z), we find that the 
Majorana bilinear along the (1,1,1) direction creates a 
pair of fundamental excitations: 

E{x,y,z) ■ =1+^. (34) 

Therefore, the fundamental excitation in this model is 
clearly not a fracton. We now demonstrate that the 
fundamental excitation may only hop freely along the 
( 1 , 1 , 1 ) direction, without creating additional excitations. 
Consider the variety V (/) defined by the stabilizer ideal 
I{S) = {1 E X ^ y + z ^ xy ^ yz ^ xz, xyz -h x + ^ + z + 
xy-\-yz-\- xz), i.e. the zero-locus of the generators of the 
ideal over an extension of F 2 . The following is a point on 
the variety: 

ix,y,z) = (^t, (35) 

with t in an extension of F 2 . As a result, if 1 + x^y^ G 
I{S), we must have from (35) that = (1 + for 
infinitely many t. This can only be true if n = m = 
0. As a result, the fundamental excitation cannot hop 
freely in the x^-plane. As the generators of the ideal are 
symmetric under exchanging any pair of variables (e.g. 
X i —^ y), we conclude that 1 + y'^z'^, 1 + x'^z'^ G I{S) 
only if n = m = 0 , so that the fundamental excitation 
cannot freely hop in the yz- or x 2 :-planes. From these 
results, we are led to the conclusion that 

1 + a;" 2 /™/ ^ I{S) (36) 

when n, m and i are distinct. Therefore, we have shown 
that the fundamental excitation in the model defined by 
f^{x,y,z) is restricted to hop along the ( 1 , 1 , 1 ) direction 
of the cubic lattice without creating additional excita¬ 
tions. 
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Appendix A: Stabilizer Excitation Maps in 
Polynomial Representations of Majorana 
Hamiltonians 


In this section, we prove a statement presented in the 
main text, that in a commuting Majorana Hamiltonian, 
of the general form 


m 


(Al) 


two Majorana operators On^ and commute over all 
lattice sites (i.e. [On\Om^] = 0 for any n, m) if and 
only if their respective stabilizer maps Si and Sj satisfy 
Si • Sj = 0. Furthermore, we show that the excitation 
map Ei = Si may be constructed for each operator in 
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the Hamiltonian, so that, given the polynomial represen¬ 
tation P of some operator, the quantity Ei • P yields a 
polynomial encoding the pattern of excitations created 
by that operator when acting on the ground-state. 

Consider a d-dimensional lattice with translation 
group A = Z^, and a single Major ana fermion per lat¬ 
tice site. At a given site, the identity operator 1 and the 
Majorana fermion 7 form the group Z 2 = Z/2Z under 
multiplication since 

lx 7 = 7 xl =7 (A2) 

lxl= 7 X 7 = l. (A3) 

We refer to this as the ‘Majorana group’ Ad at a given 
lattice site. The group A4 naturally forms a vector space 
over the finite field F 2 since for some a G Z 2 and m ^ M 
we may define the Z 2 action a • m = G M. Further¬ 
more, any element of the translation group g ^ A has a 
natural action on A4 by multiplication. For instance, let 
c,d G Z 2 , so that we may represent the action of g on 
the operator ( 7 )^ simply as • c G F 2 . A more general 
operator, for instance {g ’ jY ^ {h • jY written 

as g • c-\- h • d ^ F 2 . 

Now, we consider a d-dimensional lattice with an n-site 
basis so that each site again contains a single Majorana 
fermion; n is restricted to be an even integer so that 
we may have a well-defined number of eomplex fermions 
per lattice site. We represent each of the n Majorana 
fermions at a given lattice site as 7 ^, with j = 1 , ..., 
n. Recall that distinct Majorana fermions anti-commute. 
Let ^ G A = Z^ be an element of the lattice translation 
group. From the natural action of the translation group 
on the Majorana fermions, we may write any two Her- 
mitian operators Oi and O 2 up to overall factors of Pi 
as 


OisII 

geA 

geA 


n (S' ■ 

j=i 


Cjig) 


n (S' ■ 'yp 

j=l 


djig) 


(A4) 

(A5) 


with the coefficients Cj{g), dj{g) G Z 2 . On physical 
grounds, we require that each operator is the product 
of an even number of Majorana fermions so that the to¬ 
tal fermion parity is conserved. This is equivalent to the 
condition that 

n n 

ZZ‘^i(S') = ZZ^i(S') = 0 ^ 1^2 (A6) 

geA 3 = 1 geA 3 = 1 


Two operators commuting with the total fermion par¬ 
ity commute if and only if they have overlapping support 
on an even number of Majorana fermions. As a result, 
the operators Oi and h • O 2 commute for some h G A if 


n 

= 0 . (A7) 

geA j=l 


The condition (A7) may be simply imposed by repre¬ 
senting the operators Oi and O 2 as 


III 

M 

( Cl ( 5 ) \ 

C2(5) 

^2 = 'pg- 

( di{g) \ 
d2{g) 

(A 8 ) 

geA 

\C"(5) / 

geA 

\ dn{g) j 



respectively. We further define di by combining the an¬ 
tipode map g ^ g~^ with transposition of the Si vector. 
For example, 

S'! = Z]s'”^ • ( ^”(s') ) (^9) 

geA 

Now, we demonstrate the following statement: S 2 'Si = 0 
if and only if [Oi, i • O 2 ] = 0 V G A. 

We demonstrate this by explicit calculation. First note 
that: 


n 

S 2 -Si= y] '^h-^g-Cj{g)dj{h) 

g,heA j=l 

n 

ieA geAj=l 


(AlO) 

(All) 


We now observe that S 2 'Si = 0 if and only if the quantity 
in brackets vanishes \/ i e A. However, this is precisely 
the condition (A7) required so that [Oi, £'02] = 0. This 
completes the proof. 

A natural consequence of this proposition is the follow¬ 
ing. Consider N types of operators at each lattice site, 
and a Hamiltonian of the form: 


H = - 




(A12) 


where the sum is over all lattice sites. We may now rep¬ 
resent each operator at a particular lattice site by a 
stabilizer map S^. Then, all of the operators appearing in 
the Hamiltonian mutually eommute (i.e. [On\Om^] = 0 
V i, j, m, n) if and only if 

^i'Sj=0 (A13) 


for any pair i, j. 

From the stabilizer map, we may determine how the 
action of an arbitrary operator O on the ground-state cre¬ 
ates a pattern of excitations. Let S be the stabilizer map 
for an ideal Majorana Hamiltonian with a single operator 
per lattice site and an n-site unit cell. Furthermore, let 
P be the polynomial representation of some operator O. 
We represent S and P as 



( Cl ( 5 ) \ 
C2{g) 

73 

III 

M 

( di{g) \ 
d2{g) 

(A14) 

geA 

\ Cn{g) J 

geA 

\ dn{g) j 
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respectively. Recall that the condition (A7) determines 
whether the operator O commutes with a given stabi¬ 
lizer appearing in the Hamiltonian. Whenever, O anti¬ 
commutes with a stabilizer operator, it creates an exci¬ 


tation at the location of the stabilizer. The pattern of 
excitations created by O is then specified by the expres¬ 
sion 


n 

heA geAj=l 


h,geA j=l 


Therefore, we define the excitation map 


E = S 


(A16) 


SO that E • P yields a polynomial representation of the 
pattern of excitations created by the operator O when 
acting on the ground-state. 


Appendix B: Stabilizer Maps for Ideal Majorana 
Hamiltonians with a Two-Site Basis 


We now prove the following statement: a topologically- 
ordered, ideal Majorana Hamiltonian with a two-site ba¬ 
sis and a single operator per lattice site, is specified by a 
stabilizer map is of the form: 


/ f{xi,...,Xd) \ 

5= _ (Bl) 

\f{xi,...,Xd) / 


with f{xi ,..., Xd) G R = F 2 [xi,..., . Note that S is 

uniquely defined only up to monomial multiplication on 
each of its entries. 

We demonstrate this as follows. For an ideal Majo¬ 
rana Hamiltonian with a two-site basis, let the stabilizer 
map be of the form 5'^ = (/, ^) with /, ^ G R. Recall 
that the excitation map is given hy E = {f,g). We re¬ 
quire on an open lattice that ker(F^) = im(S'), so that 
any degenerate ground-states of the Hamiltonian are lo¬ 
cally indistinguishable. Note that ker(F^) and im(5') are 


dj{h) = 

i=i 


geA heA 


= S-P 


(A15) 


defined as 


ker(£’) = 
im(S') = <j 


a 

p 

c-f 
c-g 


af = pg 
c&r' 


(B2) 

(B3) 


Let g and / take the form g = hg' and / = h/', 
for some h G R, so that gcd(/, = h, and so that 
gcd(^', /') = !. In this case, we observe that 


9 ' 

7 


G ker(F^) 


(B4) 


If this element is generated by im(5'), it must be the case 
that 



for some c e R. However, this equation is only satisfied 
\i c-h-c - h = 1, which implies that both c and h must be 
monomials. As S is only uniquely defined up to mono¬ 
mial multiplication on its entries, we may restrict the 
remainder of our analysis to stabilizer maps S = (/, g7 
with / and g co-prime so that gcd(/, = 1. 

Since / and g are co-prime, the condition af = jSg 
for a vector (a,/3)^ to be in ker(F^) is satisfied only if 
a = c-g and /d = c• / for some c e R. Now consider that 
in order for im(5') = ker(F^) on an open surface, we must 
have that 



for some d e R. This is indeed satisfied if c = d, so 
that f = g. Now, we have that the most general form 
of the stabilizer map (up to monomial multiplication on 
each entry) is S = (/,/)^ for some f e R such that 
gcd(/, /) = 1. The commutativity condition S ' S = 0 is 
trivially satisfied. This completes the proof. 











